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. ABSTRACT 



In September 2011, the Herschel Space Observatory performed an observation campaign with the PACS photometer observing the 
asteroid (101955) 1999 RQ36 in the far infrared. The Herschel observations were analysed, together with ESO VLT-VISIR and 
^ . Spitzer- IRS data, by means of a thermophysical model in order to derive the physical properties of 1999 RQse- We find the asteroid 

■ has an effective diameter in the range 480 to 51 1 m, a slightly elongated shape with a semi-major axis ratio of a/b=1.04, a geometric 
^ ■ albedo of 0.045^[5 jj^, and a retrograde rotation with a spin vector between -70 and -90° ecliptic latitude. The thermal emission at 

, wavelengths below 12/jm -originating in the hot sub-solar region- shows that there may be large variations in roughness on the 

surface along the equatorial zone of 1999 RQ36, but further measurements are required for final proof. We determine that the asteroid 
has a disk-averaged thermal inertia of r=650 Jm^-s^^-^K"' with a 3-cr confidence range of 350 to 950 Jm^-s^^-^K"', equivalent to what 
is observed for 25143 Itokawa and suggestive that 1999 RQ36 has a similar surface texture and may also be a rubble-pile in nature. 
^ , The low albedo indicates that 1999 RQ36 very likely contains primitive volatile-rich material, consistent with its spectral type, and 

■ that it is an ideal target for the OSIRIS-REx sample return mission. 

!>■■ 

^■f^ . Key words. Minor planets, asteroids: individual - Radiation mechanisms: Thermal - Techniques: photometric - Infrared: planetary 

1^-^ ■ systems 

d- 

^ '. 

' 1 . Introduction belt and are believed to be primitive and volatile-rich. However, 

,—1 ■ the Polana family in the inner belt contains B-types, and dynam- 

. On May 25, 2011, NASA announced the selection of OSIRIS- j^^j ^^^^-^^ ^^^^ ^^^^^ rq^^ ^- ^^ ^ Hberated 

_ > ■ REx as the third selected mission of its New Frontiers Program. ^^^^^^ ^^^-^ ^^^-^ considered as an important inner-belt 

^ . The OSIRIS-REx mission will be launched in September 2016, ^^^^^^ low-albedo NEAs (Campins et al. M Cellino et al. 

^ . fly by the Earth for a gravity assist in September 2017, and 200I ). Additionally, since its orbit makes it especially accessible 

; encounter the neai-Earth asteroid (NEA) (101955) 1999 RQ36 ^o spacecraft and it has been identified as a potentiafly hazai'dous 

■ - ■ (Lauretta et al. EHISll m November 2019. Proximity operations ^^^^^^-^ ^j^^^^- 2009 ), it was considered an ideal choice 

at the asteroid last flirough Mai'ch 202 1 , when it will acquire up ^^3^,^ OSIRIS-REx sample return mission (Lauretta et al. 

to 2 kg sample of its surface material. The sample return capsule 2012 ) 
will return to the Earth's surface in September 2023. 

Based on its visible spectrum, 1999 RQ36 is classified as a R^^^^" ^^8^^' ^^^en during its previous two closest ap- 

B-type and CM chondrite meteorites are considered to be the f^^'^hes in 999 and in 2005 (and the next occiirred in 

corresponding spectral analog (Clark et al. 201 1 ). B-class aster- ^f^f}^"^^^' ^01 1 when our observations were performed) found 

oids are found mostly in the middle and outer regions of the main 1^99 RQ36 tobean irregular spheroid about 580 m in diameter 

(Nolan et al. 120071. The nearly spheroidal shape suggests that it 

* Herschel is an ESA space observatory with science instruments might be a strengthless rubble pile, formed by mechanisms sim- 

provided by European-led Principal Investigator consortia and with im- ilar to near-Eaith binary systems (Walsh et al. 2008 ). The latest 

portant participation from NASA. analysis, where all available radar data was combined with vis- 

** Analysis is also based on observations collected at the European ible hghtcurves (Nolan et al. I2012I I. derived a mean diameter 

Southern Observatory, Chile; ESO, No. 287.C-5045(A) of 493 + 20 m (mean equatorial diameter of 545 + 15 m) and a 
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"spinning top" shape similar to 1999 KW4 (Ostro et al. I2006I I. 
but with a less well-defined equatorial ridge. 

The objec t's rot ation period is 4.297812 + 0.000001 hours 
(Nolan et al. I2012I I. which was derived from combined radar 
and lightcurve measurements. This period is very similar to the 
4.288 h found by Nolan et al. (j2007 ) from radar measurements 
alone and to the 4.2968 + 0.0017 h based on the analysis of 
lightcurves (Hergenrother et al. 120121) and about twice the value 
given in Krugly et al. (2002) derived from incomplete R-band 
lightcurve observations during the 1999 opposition. The rota- 
tion period suggests that it has not been greatly spun up by tidal 
or radiation forces. However, the radar-derived pole orientation 
{Ppgi^ - -90° + 15°) is the e quilibrium state for the YORP effect 
(Vokrouhlicky et al. 120031 ). The negative pole orientation indi- 
cates a retrograde rotation. 

Emery et al. (2010) find a size of 610m and estimate its ther- 
mal inertia to be around 600 Jm"^s"*''^K"', assuming a moderate 
surface roughness (/:m..s'.-slopes ~20°). Their results are based 
on Spitzer-IRS spectra and a dual-band thermal lightcurve, taken 
on May 3/4, 2007 at a phase angle of about 63°. 

In this paper, we first present our R4CS and ESO-VISIR 
Director Discretionary awarded Time (DDT) observations taken 
of 1999 RQ36 and their data reduction. We follow this with 
the description of the Spitzer-IRS measurements and the corre- 
sponding data reduction. We analyse the derived flux densities 
by means of our thermal model. We proceed to present our ra- 
diometrically derived properties (Section[3]). In Section|4]we dis- 
cuss the object's size, shape, and albedo in a wider context and 
study the influences of surface roughness and thermal inertia in 
detail. We finally conclude the paper with a summary of all de- 
rived properties and with the implications of our results for the 
preparations for the OSIRIS-REx mission. 



2. Observations of 1999 RQ36 

2.1. Herschel-PACS observations and data processing 

The European Space Agency's (ESA) Herschel Space 
Observatory (Pilbratt et al. 120101 ). launched in 2009, per- 
forms observations from the 2'"' Lagrangian point (L2) at 
1.5 10^ km from Earth. It has three science instruments on 
board covering the far-infrared part of the spectrum, of which 
the PACS Photometer (covering the wavelength range of 
60-210jum) has been used to observe 1999 RQ36. 

The PACS photometer (Poglitsch et al. 120101 ) imaged the 
asteroid on September 9, 2011 (between 19:00 and 21:00 UT) 
at a heliocentric distance of 1.0146-1. 0144 AU and a Herschel- 
centric distance of 0.1742AU, and the phase angle was 85.3- 
85.4°. The time used for the observation was DDT-awarded 
specifically due to this optimum closest approach to Earth. The 
70/160/im scan-map observation was repeated ten times (5 times 
at scanangle 70° and 5 times at 110°), the 100/160/im observa- 
tion was repeated 14 times (7 times in each scanangle). 

The PACS scan-map measurements were processed using 
HIPE7.0 (Ott et al. 120091 ). The reduction steps were very sim- 
ilar to those used for trans-Neptunian objects (Santos-Sanz et 
al. (2012; Mommert et al. 2012 Vilenius et al. l20T2b . includ- 
ing corrections for moving targets. The final, single-repetition 
PACS maps (10 at 70//m, 14 at 100 //m, and 24 at 160 yum) have 
been median-averaged to remove (or at least decrease) the other- 
wise large instrumental noise. The background confusion noise 
has been eliminated as much as possible via techniques that 
have been developed for other moving sources (see detailed de- 
scription in the electronic appendix in Santos-Sanz et al. 120121 ). 



The derived fluxes were aperture- and colour-corrected to ob- 
tain monochromatic flux densities at the PACS reference wave- 
lengths. The colour correction values for 1999 RQ36 of 1.01, 
1.02, 1.07 in blue (70//m), green (100 //m) and red (160;um) 
bands are based on a thermophysical model SED (spectral en- 
ergy distribution), corresponding roughly to a 250 K black-body 
curve (PACS photometer passbands and colour correction fac- 
tors, 201 1 ). The flux calibration was verified by a set of five high- 
quality fiducial stars (ySAnd, aCet, aTau, a Boo, and yDra), 
which were observed multiple times in the same PACS observing 
mode (PACS photometer - point source flux calibration. 1201 II ). 

The final values obtained for the PACS instrument datasets 
were deemed to be accurate within 5% based upon existing cal- 
ibrations. The absolute flux calibration uncertainty is included 
in the values given in Table [1] (added quadratically to obser- 
vational errors and errors related to aperture photometry). The 
fluxes given in Table[T]are colour-corrected, monochromatic flux 
densities at 70.0, 100.0, and 160.0/im, the three PACS reference 
wavelengths for the blue, green, and red bands, respectively. 
Figure[T]shows the combined PACS images aligned with the out- 
put from our thermal model of what PACS has observed during 
its observation campaign. 




PACS view: 09-Sep-2011 19:55 UT; T [K] 




300 



250 



200 



150 



Fig.l. The combined PACS images (top), at 70 (blue), 100 
(green), and 160yum (red), along with the output from our ther- 
mal model (bottom) of 1999 RQ36, at the geometry of Herschel's 
viewpoint. The shorter wavelengths effectively sampled the 
warmer facets of 1999 RQ36's surface, located in the subsolar 
region. Note that Herschel cannot spatially resolve the target. 



2.2. ESO-VISIR observations and data processing 

Besides the DDT awarded for Herschel, we also were awarded 
DDT to observe 1999 RQ36 in September 2011 via ground- 
based N- and Q-band observations with the ESO-VISIR instru- 
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Table 1. Observation summary. 



Start Time UT 


Duration 


Julian Date 




FD 


FD,,, 




^obs 


a 


Observatory/ 


yyyy-mm-dd hh:mm:ss 


[sec] 


mid-time 


[pan] 


[mJy] 


[mJy] 


r [AU] 


[AU] 


[deg] 


Instrument 


2007-05-03 00:00:14 


Start 


2454223.50016 








1.1237961 


0.5056454 


-63.5 


Spitzer-PUI 




llx 254 


2454224.23335 


16/22 


Lightcurve 


1.1262212 


0.5089081 


-63.3 


Spitzer-PUI 


2007-05-04 11:15:50 


End 


2454224.96655 








1.1286398 


0.5122203 


-63.1 


Spitzer-PUI 


2007-05-03 21:39:13 


4261 


2454224.42689 


5.2-38 


SED 


1.1268605 


0.5097779 


-63.2 


Spitzer-IRS 


2007-05-03 23:45:57 


4260 


2454224.51489 


5.2-38 


SED 


1.1271494 


0.5101723 


-63.2 


Spitzer-IRS 


2011-09-09 19:01:12 


1414 


2455814.30068 


70.0 


27.0 


1.7 


1.0146021 


0.1742221 


-1-85.3 


Herschel/PACS 


2011-09-09 19:59:50 


1414 


2455814.34140 


70.0 


24.3 


1.5 


1.0144716 


0.1742140 


+85.4 


Herschel/PACS 


2011-09-09 19:25:49 


1978 


2455814.32104 


100.0 


14.2 


1.1 


1.0145368 


0.1742180 


-^85.4 


Herschel/PACS 


2011-09-09 20:24:27 


1978 


2455814.36176 


100.0 


12.0 


1.2 


1.0144064 


0.1742101 


+85.4 


Herschel/PACS 


2011-09-09 19:25:49 


1978 


2455814.32104 


160.0 


6.5 


2.3 


1.0145368 


0.1742180 


+85.4 


Herschel/PACS 


2011-09-09 19:59:50 


1414 


2455814.34140 


160.0 


4.4 


1.9 


1.0144716 


0.1742140 


+85.4 


Herschel/PACS 


2011-09-09 20:24:27 


1978 


2455814.36176 


160.0 


7.6 


2.0 


1.0144064 


0.1742101 


+85.4 


Herschel/PACS 


2011-09-17 09:21:59 


4x ~2min 


2455821.89306 


8.59 


23.7 


3.6 


0.9908803 


0.1798955 


+89.4 


ESO VLT/VISIR 



Notes. The eleven Spitzer-IRS Peak-up Images (PUI) have the AORKEYs 2141 2864/4400/4912/3120/5168/8240/3632/4144/4656/3888/3376. 
The data contain simulataneous 16 and 22yum measurements. The re-reduced Spitzer-IRS data have the AORKEYs 21412352 and 21412608. 
They are separated by about half a rotation period. See also Emery et al. ( 2010J for more observational details and the presentation of preliminary 
results. The Herschel-PACS observations have the observation identifiers 1342228379 ... 1342228382, taken on Operational Day 849, the ESO 
VLT-VISIR observations are related to the programme-ID 287.C-5045(A). 



ment (Lagage et al. 120041) mounted on the 8.2 m VLT telescope 
MELIPAL (UT 3) on Paranal. 

The service-mode observers attempted twice to execute our 
four observing blocks in imaging mode, including each time the 
PAHl (8.59yum), Nell (12.81jum), and the Q2 (18.72 ^um) fil- 
ters. The 1999 RQ36 observations were done in parallel nod- 
chop mode with a throw of 8". Owing to poor observing con- 
ditions and difficulties during the target acquisition phase, com- 
bined with end-of-night constraints, the science observations 
were never executed. Nevertheless, we found 1999 RQ36, with 
both the positive and two negative nod-chop beams on-chip, in 
the PAHl target acquisition images taken on September 17, 20 11 
at an airmass of 1.79. We processed the four useful raw images 
using the ESO VISIR pipeline Gasgan^ to produce a combined 
image suitable for aperture photometry. It was possible to cali- 
brate those data using the standard star HD 217902 (9.923 Jy at 
8.59/im, taken at airmass 1.30) that was observed before the ac- 
quisition of 1 999 RQ36 and another one, HD 1 68 1 5 ( 1 1 .679 Jy at 
8.59yL(m, taken at airmass 1.40), which was observed the day af- 
ter in the same filter and observing mode. Based on the derived 
count-to-Jy conversion factors and the experience from large cal- 
ibration programmes (Schiitz & Sterzik I2005I I. we were also 
able to coiTect for the airmass difference between the calibra- 
tion stars and our science target. The counts decrease by about 
10% for an airmass increase of 0.5. The fluxes of the standard 
stars and our target were extracted in an identical way to stan- 
dard aperture photometry techniques. In the case of 1999 RQ36 
we checked that the first nod-chop pair was consistent with the 
second (i.e. that all three beams were in all the acquisition im- 
ages) by also reducing them separately. The derived fluxes from 
each pair were consistent with those produced by combining 
both pairs together. We applied a colour correction of 1.004 
(for 1 50-250 K blackbodies, considering the full atmosphere and 
bandpass transmission profiles) to the measured 1999 RQ36 flux 
to obtain a mono-chromatic flux density at the PAHl reference 
wavelength (see Table[T]l- For the error calculation (E. Pantin, M. 
van den Ancker & R. Siebenmorgen, priv. comm.), we quadrati- 



caUy added the following error sources: error in count-to-Jy con- 
version factor (2%), error of stellar model (2%), photometry er- 
ror (10%), error of colour correction (2%), estimated additional 
uncertainty for measurements at airmass larger than 1.5 (5%), 
and flat-field error for an off-centre position in the acquisition 
images (10%), summing up to a total of 15% eiTor in the derived 
flux density. 



2.3. Observations of 1999 RQ36 with fhe Spitzer Space 
Teiescope 

To complete the dataset, we also analysed and utilised data from 
the Spitzer Space Telescope (Werner et al. '2004V Emery et al. 
(2010) present part of these observations, but the results were 
based on a preliminary data reduction alone. 

Spitzer observed 1999 RQ36 in different modes. Thermal 
IR spectra from 5.2 to 38;um were taken with the Infrared 
Spectrograph (IRS) of opposite hemispheres of the body 
on 3/4 May 2007 21:42:03 - 22:46:57 (ID 21412352) and 
23:48:51 - 00:53:40 UT (ID 21412608). Photometry at 16 
and 22yum was obtained with the IRS peak-up imaging (PUI) 
mode (IDs 21412864/4400/4912/3120/5168/8240/3632/4144/- 
4656/3888/3376), taken at regular intervals between 3 May 2007 
00:00:15 UT and 4 May 2007 11:11:37 UT The 11 measure- 
ments cover the fuU rotation period uniformly. The phase an- 
gles, the heliocentric, and Spitzer-centric distances are given in 
TablelU 

The Spitzer-IRS observations were reduced with SPIC^ in 
the recommended way following the data analysis cookbook 
for moving target^ We used the PBCD-productfl from the 
Spitzer Heritage Archive SHA0, subtracted the background mea- 



' |http : //www ■ eso ■ org/ gasgzino | 



- http : //i rsa . ipac . caltech . edu/data/SPITZER/docs/dataanalysi 
[May 7"', 2012] 

' Ihttp : //irsa . ipac . caltech . edu/dat a/SPITZER/docs/dataanalysi 
[May 7*, 2012] 

Post-Basic Calibrated Data: Level 2 data 
' jhttp : // sha . ipac . caltech . edu/applications/Spitzer/SHA/ j 



[May 7^2012] 
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surements from the on-source measurements, and produced two 
spectra for each order in the short-low (SL) and long-low (LL) 
module^ using the appropriate uncertainty and mask files. We 
averaged the resulting A-B and B-A spectra. The low flux in the 
SL2 range (5.2 to 7.5 fim) was averaged using a boxcar filter; in 
addition, we resampled the data in bins of uniform resolution R 
= A/AA, with R = 40 (SLl: 7.6 to 14.3 A^m) and R = 30 (SL2) 
to better measure the continuum signal of 1999 RQ36 and to re- 
duce the number of data points for the thermophysical model 
(TPM) analysis. For each bin we determined the weighted av- 
eraged wavelength, the median flux, and the root-sum-square of 
the propagated errors, including a 5% absolute flux calibration 
error as specified by the Spitzer/IRS team. Our derived fluxes 
and errors agreed within the given errorbars with the published 
figures in Emery et al. ( 1201 Oi l and are shown in Fig.|2] 




Wavelength [urn] 

Fig. 2. Both Spitzer IRS spectra (SL-part only) are shown to- 
gether with 2 PUI measurements taken close in time at similar 
rotational phases. Our best TPM solution is shown as a solid line. 



2.4. Observation summary 

In total, including the Herschel, ESO-VISIR and Spitzer dataset 
(see Table[TJ, we were able to place into our model a set of obser- 
vations covering the wavelengths from 6 to 160 fim, phase angles 
of -63°, -1-85°, and -h89° (positive before opposition, negative af- 
ter opposition), and covering nicely all rotational phases. For the 
TPM analysis (Sect. O we use the true observing and illumina- 
tion geometries (Spitzer-, Herschel-, VLT-centric) for each mea- 
surement individually as given in Table[T] The importance of this 
dataset becomes clear in the coming section where it can be seen 
that our comprehensive coverage of wavelength and phase angle 
is key to extracting the main thermal parameters output from our 
thermal model. 



* Note, that we used only the SL-part of the observations. The LL 
part is more noisy and does not constrain our analysis in a significant 
way. 



3. Thermophysical analysis of the observational 
data 

We started the thermophysical model (TPM) analysis of the 
complete dataset by using a simple spherical shape model. For 
the object's rotatio n period we used 4.297812 + 0.000001 hours 
(Nolan et al. l2012l i. Thermal data are very sensitive to the orien- 
tation of the spin axis (see e.g., Muller et al. 201 1 ) and we used it 
as a free parameter in our analysis. Due to computing speed lim- 
itations we restricted this exercise to 60 different orientations for 
the rotation axis distributed over the entire sphere. The ecliptic 
latitude and longitude pairs to describe the spin vector orienta- 
tion correspond to the 60 vertices of a mathematical description 
of a truncated icosahedrorQwith 12 regular pentagonal faces and 
20 regular hexagonal faces. To complete the picture, we added 
a few key spin axis orientations: pro- and retrograde rotations 
with the spin vector perpendicular to the ecliptic plane (JSpoie = 
+90°), pole-on viewing geometries as seen during the Spitzer, 
the Herschel, and the VISIR observations, and another four ori- 
entations (with /3p„ie < -70°) to explore the YORP spin-up equi- 
librium state (Vokrouhlicky et al. I2003I I. which corresponds to 
the latest value derived by Nolan et al. (2012): The obliquity 
180° case for the spin vector (J3p„ie = -83.98°, Apoi, = 92.17°) 
is very well covered by by this set of different spin axis orien- 
tations, which included values of -78°, -80°, and -90° for jSpoie 
(each time for a range of Apoie values). 

The spherical shape model with the almost 70 different spin 
axis orientations were imported into our TPM code (Lagerros 
119961 119971 119981 1. Details about the tec hniques have recently 
been described in Muller et al. (201 1) and O'Rourke et al. 
(2012). The TPM produces accurate thermal IR spectra and ther- 
mal lightcurves, taking a number of physical and thermal pro- 
cesses into account. In the TPM, the object is described by a 
given size, shape, spin-state, and albedo, placed at the true ob- 
serving and illumination geometry. The TPM considers a 1-d 
heat conduction into the surface and allows surface roughness to 
be included, described by "f , the fraction of the surface covered 
by spherical crater segments and "p", the r.m.s. of the surface 
slopes, connected to the crater width-to-depth ratio (Lagerros 
ITW8 ). A default setting (Muller et al. |T999] l wifli p=0.7 and 
f=0.6 has been used as baseline. 

The contributions of the subsurface emission at longer wave- 
length is accounted for by a wavelength-dependent emissivity 
decreasing from 0.9 at mid/far-IR (5 to 100 ;um) to about 0.8 
in the submillimetre/millimetre range of the spectrum, derived 
from a combined set of large main-belt asteroids (Muller and 
Lagerros ll9981l2002l l. Here, all radiometrically critical measure- 
ments are taken at wavelength below 100 yum, and there is no dif- 
ference in using a constant emissivity of 0.9 or the wavelength- 
dependent emissivity model presented in Muller & Lagerros 
(l2502l l. 

For the full treatment of the energy balance on each sur- 
face facet within the model, we need to describe the amount 
of reflected light, described by the commonly used H-G-system 
(Bowell et al. I19891 I. The value for the absolute V-band mag- 
nitude Hv was determined by Krugly et al. (120021) and re- 
cently refined by Hergenrother et al. (120121 1. The best value of 
Hy = 20.65 mag was obtained by a linear phase function fit 
(phase slope /?= 0.039+0.005 mag) to V-band photometry cover- 
ing the phase angle range from 15° to 100°. Hergenrother et al. 
('2012 ) discuss the issue of a shallow opposition effect for low- 
albedo asteroids and conclude that a more realistic Hy-value for 



' |http : //en ■ wikipedia . org/wiki/Truncated-icoszihedron| 
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1999 RQ36 would be 20.5+0.3 mag. For G (slope parameter) we 
used the Ki-ugly et al. ("2002^ value of 0. 12. 

Based on all available thermal data we searched for the best 
size and albedo solutions via well-established radiometric tech- 
niques (e.g., Harris & Lagerros 2002 and references therein) by 
using the above described settings. For each of the above spin 
vector orientations we calculated the reduced ^^-value for a wide 
range of typical thermal inertias (Delbo et al. 12007 1. We calcu- 
lated the the reduced x^-^^lnes in the following way: xleduced ~ 
l/(N - 2) * Y^iiflxobs - .flx„,„d)/errohs)^, with N the number of 
degrees of freedom. Figure [3] shows the result of the analysis. 
The lowest values are found for retrograde rotations (/3p„ie < 
-70°) and thermal inertias around 600 Jm^^s "^K The ther- 
mal inertia 3-cr confidence interval is from 400 to just above 
lOOOIm^^s^^-'^K"'. The 3-cr confidence interval 2 < + 
is based on actual - {{obs - mod)/err)^ calculations (not 
the reduced as shown in Fig. |3]l. 



-i 1 1 — 1 — 1 — i i i i I 1 1 — 1 — 1 — I I I i 

- spherical shape model 




: pv= 0.033-0.060 



j 1 1 1 1 ■ ■ ■ ■ I 1 1 1 1 ■ ■ ■ ■ I 1 I 

10 100 1000 

Thermal inertia [J m"^ s"°'' K"'] 

Fig. 3. t6st for the full set of spin vector orientations and 
the full range of thermal inertias. The retrograde solutions with 
l3poie < -45° are given as solid lines, the prograde ones with 
l3poie > +45° as dashed lines, the ones where the latitude of the 
rotation pole is in the range -45° < fipoie < +45° are shown 
in dotted lines. The thicker solid and dashed lines represent the 
/3poie = ±90° solutions, the thick dashed-dotted line the Herschel 
pole-on geometry. 

Considering this possible F-range in combination with a 
spherical shape model and a pole axis pointing to the ecliptic 
south pole (j3poie = -90°, retrograde rotation) leads to a size in 
the range 498 to 513 m and a geometric albedo pv in the range 
0.043 to 0.046. The true eiTors are larger: (i) the TPM fit to the 
observations is not perfect (mainly due to the spherical shape 
model), and the r.m.s. of all observation/model ratios (see Figs.|4] 
& |5]l has to be included (dominating error contribution for the 
diameter); (ii) the solution also depends slightly on the selected 
spin vector orientation within the possible range /3poie < -70°; 
(iii) the H-magnitude is given with an error of +0.3 mag (dom- 
inating error contribution for the albedo). Combining all errors 
quadratically led to an effective size of 503^}4m for a spherical 
shape and a geometric albedo of 0.045;^q q}2. These numbers are 
based on the assumption that 1999 RQ36 has a "default rough- 
ness" that is typical of large main-belt asteroids (see also discus- 
sion in Section|4]below). 



4. Discussions 

4.1. Spin-axis orientation, sense of rotation 

The before/after opposition data are the key element for extract- 
ing the sense of rotation. Our Herschel observations are all taken 
before opposition, i.e. leading the Sun (here positive phase an- 
gles). This means that Herschel has observed the target, while 
the unilluminated part of the surface was cold (see Fig. [1]) cor- 
responding to the pre-dawn part of the surface (for a retrograde 
sense of rotation). Spitzer has seen 1999 RQ36 after opposition 
when the unilluminated part was still warm after it had rotated 
out of the Sun (see Fig.|7]i. The importance of the sense of rota- 
tion can clearly be seen in Fig. [3] The prograde rotation options 
produced very high values for the reduced x^ and can therefore 
be excluded with high confidence. The models with spin axes 
close to the ecliptic plane (including pole-on geometries) con- 
strain the thermal inertia only very little; nevertheless, the re- 
duced -values are very high so these solutions can be ruled 
out. The TPM analysis clearly favours an axis orientation with 
/3poie < -70° in the ecliptic coordinate system. This agrees with 
the findings by Nolan et al. (2(312) that the pole is close to the 
equilibrium state for the YORP effect (Vokrouhlicky et al. l2003l l. 
1999 RQ36 has an orbit inclination of 6.02°, and Q, the argu- 
ment of ascending node, is 2.17°. The YORP predicted asymp- 
totic zone with obliquity e of 180° would therefore correspond to 
/3poie - -83.98° and Apo/g = 92.17° in our convention, very close 
to the south ecliptic pole. 
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Fig. 4. Ratio between the observed fluxes and the correspond- 
ing TPM predictions as a function of phase angle (negative: af- 
ter opposition; positive: before opposition). The model calcu- 
lations are based on a thermal inertia of 600Jm"^s " ''K^' and 
p-0.2 for the surface roughness. The VISIR data point can only 
be matched by assuming a higher surface roughness (see also 
Section l4~2l for further discussions). 



4.2. Surface rougliness 

Due to the lack of observations close to opposition, our dataset 
does not constrain the surface roughness that plays a more dom- 
inant role at phase angles close to opposition very well (Muller 
2002). In Section [3] we used the "default roughness" that was 
derived for large main-belt asteroids that are covered with a low 
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thermal conductivity regolith (Muller et al. I1999I I. To study the 
influence of roughness on the final size and albedo solutions, 
we analysed all thermal data, but this time for diff'erent levels of 
roughness. For simplicity, we assumed that 100% of the surface 
is covered by craters (f=1.0), which is a reasonable assumption 
considering that the model crater definition also includes micro- 
craters down to very small scales where the geometric optics 
approximation is still valid (Muller & Lagerros 1998). We var- 
ied the r.m.s. of the surface slopes p in steps of 0.1 from a per- 
fectly smooth surface (p=0.0) to a very rough surface {p-\.0) 
and repeated our radiometric analysis. The reduced ;t'^-minima 
indeed reach acceptable values for almost all levels of rough- 
ness: A very low surface roughness with p < 0.1 would lead to 
a thermal inertia of 450Jm"^s "^K ' which corresponds to an 
effective size of 490m and a geometric albedo of p^ = 0.047. A 
surface roughness with p > 0.9 would point to a thermal inertia 
of 850Jm-2s-0-5K-i with D,// = 508 m and pv = 0.044. This 
size range is entirely consistent with the possible range given 
by radar (Nolan et al. l2012l l. but it confirms that our dataset suf- 
fers from the degeneracy between surface roughness and thermal 
inertia. Figures |4] and |5] show the comparison between obser- 
vations and model predictions assuming a roughness of p-Q.2 
(with f=1.0), similar to the moderate surface roughness (r.m.s. 
slopes -20°) used by Emery et al. (|20T0] | to fit the IRS data. 

The only data that do contain some information about rough- 
ness are the data at short wavelengths below 12yum (IRS- 
spectrum and VISIR data point). The Wien part of the SED is 
strongly influenced by the hottest temperatures in the sub-solar 
region, while data at longer wavelengths beyond the thermal 
emission peak are more connected to the disk-averaged temper- 
atures. This can be seen in Fig. |6] the low-roughness case fits 
the IRS-data very well (top of Fig. |6]l, but the model prediction 
for the VISIR data-point is overestimated by about 30%. If we 
use a high-roughness solution, then it is the other way round, 
where the VISIR data-point is well matched by the model predic- 
tion, but the IRS-spectrum divided by the corresponding model 
prediction shows a slight, wavelength-dependent mismatch (bot- 
tom of Fig. |6]l. The effects are small, and the VISIR data point 
has been taken in marginal conditions and as an acquisition 
frame. Also, modelling aspects like the implementation of sur- 
face roughness might play a role. The VISIR and IRS dataset 
are taken at very large and very different phase angles where the 
roughness effects on the IR beaming have not been studied very 
well. Nevertheless, the VISIR data point might indicate that we 
see variations in surface properties while the object is rotating. 
Based on the above specified rotation period, the two IRS spec- 
tra and the VISIR data are indeed taken at different rotational 
phases and correspond to different regions on the surface. A look 
at the reduced ;t'^-values seems to confirm this idea. The low- 
roughness model (top of Fig.|6]l matches our thermal dataset very 
well (reduced ;t'^-values below 0.75), while the high roughness 
model (bottom of Fig.|6]l clearly has problems with many of our 
data points (reduced ;^'^-values above 1.25), although it produces 
a perfect fit to the VISIR data point, which was taken at a unique 
rotational phase. But this aspect of surface roughness variations 
is clearly on vague grounds and should be confirmed via radar 
measurements or short-wavelength thermal lightcurves. 

We should note here that the thermal emission at wave- 
lengths shorter than 12/im originates mainly in the hottest sub- 
solar region, which is equatorial in the case of 1999 RQ36. Our 
speculation about possible surface roughness heterogeneity is 
therefore only related to a broad equatorial zone. Our dataset 
is not sensitive to roughness variations in the polar regions. 
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Fig. 5. Ratio between the observed fluxes and the corresponding 
TPM predictions as a function of wavelength. 



Roughness can in principle be present on any scale. Radar 
observations from Goldstone and Arecibo in 1999 and 2005 re- 
vealed a featureless surface down to the radar resolution limit of 
7.5 m (Hudson et al. 2000; Nolan et al. '2007), but so far the radar 
team(s) have not investigated whether there are also variations in 
surface roughness seen in the radar echoes (Mike Nolan, private 
communication). 

4.3. Thermal inertia F 

Pre- and post-opposition observations are the key to constrain- 
ing the thermal inertia F, even more when the observations are 
taken at large phase angles (see Fig.|4|i. But, as described above, 
uncertainty remains because of the degeneracy between surface 
roughness effects and thermal inertia. From the exercise with dif- 
ferent levels of roughness we find thermal inertias ranging from 
450 to 850 Jm"^s""-^K"'. If we also include variations for the 
spin pole (fipoie = -90° + 20°), the possible confidence range lies 
between 350 to 950 Jm"^s"''-''K"', with the most likely value at 
650Jm-2s-0-5K-'. 

A closer inspection of possible thermal inertias confirms 
that thermal inertias outside this range cause a very obvious 
mismatch between the data taken before opposition (positive 
phase angles in our convention) and the ones taken after opposi- 
tion(negative phase angles): A figure similar to Fig. |4] but with 
thermal inertias less than 350 or greater than 950Jm"^s"°^K"' 
would show a large imbalance between data taken before oppo- 
sition and data taken after opposition. The determined F-range 
is close to the 600 Jm^^s^^^K"' derived by Emery et al. (2010)1 
for 1999 RQ36 and also close to to F = 750 J m~^s~"-^K~' de- 
rived for (25143) Itokawa (Miiller et al. l2005T l. Such a thermal 
inertia value indicates regions with high conductivity, i.e. co- 
herent boulders. It can therefore be expected that the surface 
of 1999 RQ36 might show a similar texture to Itokawa's: with 
zones with different surface rock size distributions, though pos- 
sibly with fewer large boulders that would have shown up as sur- 
face features in the radar measurements (Nolan et al. l2012l l. But, 
as discussed before, our results are only connected to the warm 
region where the measured thermal flux originates. Other parts 
of the surface -at very large and very low latitudes- could have 
very different thermal properties. Seismic vibrations or YORP- 
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Fig. 6. Ratio between the observed fluxes and the coiTespond- 
ing TPM predictions as a function of wavelength. Top: assum- 
ing moderate surface roughness with a r.m.s. -slope of 0.2 gives a 
nice match to the 2 Spitzer-IRS spectra taken at opposite hemi- 
spheres (reduced < 0.75). Bottom: assuming a high surface 
roughness with a r.m.s. -slope > 0.8 gives a nice match to the 
VLT-VISIR data point, but there remains a slight trend with 
wavelengths for the IRS data which are taken at different rotation 
phases (reduced > 1.25). 



introduced regolith turn-over (see Binzel et al. '2010! and refer- 
ences therein) could have influenced these "invisible" parts by 
moving material towards the equatorial region that we "see" in 
our thermal dataset. This means that our analysis does not pre- 
clude some regions with small grains. 

4.4. Spheroidal or ellipsoidal shape 



The visual lightcurves presented by Hergenrother et al. (I2012I I 
were taken at large phase angles between 60° and 70° and 
showed a low amplitude of 0.17 magnitudes and a trimodal 
shape. It is not possible to reproduce the observed lightcurve 
shape with a simple ellipsoidal shape model of uniform albedo, 
but small shape deviations in combination with high phase an- 
gles could explain the trimodal lightcurve. The lightcurve ampli- 
tude (at a phase angle of 65°) can be reproduced by a slightly de- 
formed body with an axis ratio of a/b=1.04 and rotating around 



a c-axis that points at the ecliptic south pole. We tested if such 
a rotating ellipsoid in combination with our derived mean ther- 
mal inertia of 600 Jm^^s"*' ''K^' could also explain the two-band 
Spitzer-PUI thermal lightcurve (see also Figure 2 presented in 
Emery et al. l2010l l. Axis ratios a/b above about 1 . 1 (for the given 
spin vector and thermal inertia) can be excluded on basis of these 
16 and 22 fim lightcurves. 



16 



12 



Spitzer-PUI (Emery et al. 2010) 



■ O 16.0 /jm (blue) 
□ 22.0 /»m (red) ' 




■ TPM (a/b=1.04, r=600 SI, low roughnessp=0.1) 
I I 



100 200 

Rotational phase [deg] 

Spitzer view: 03-May-2007 00:05 UT; T [K] 




300 



I 



I 



300 



250 



200 



150 



Fig. 7. Top: Spitzer-PUI data with the TPM prediction overplot- 
ted. Bottom: The corresponding output from our thermal model 
of 1999 RQ36 as seen by Spitzer during its observation cam- 
paign. Spitzer cannot resolve the target spatially. 

Figure Q shows the result as a function of rotational phase. 
The zero rotational phase was determined in 10° steps via a 
X^-test against the observational Spitzer-PUI data. The thermal 
lightcurve of a slightly elongated shape model with a 180° obliq- 
uity of the spin vector (combined with the derived low surface 
roughness and a thermal inertia of 600 Jm^^s "^K ') explains 
the different absolute flux levels at 16 and 22yum very well and, 
at the same time, matches the two-band PUI thermal lightcurve 
and amplitude better than a simple spherical shape model. A 
lower thermal inertia would (for the given PUI observing ge- 
ometry) increase the thermal lightcurve amplitude, a larger iner- 
tia would lower the amplitude. But the uncertainties of the indi- 
vidual PUI measurements and the remaining degeneracy with 
roughness do not allow constraining the thermal inertia even 



7 



Muller et al.: Thermal & Shape properties of OSIRIS-REx Asteroid (101955) 1999 RQ36 



more. Unfortunately, the application of our PUI-derived rotating 
ellipsoidal model directly to the PACS dataset was not success- 
ful. The PACS data cover about 1 .5 h of the rotation period, but 
the quality of the PACS dataset with respect to S/N and time res- 
olution is not sufficient to repeat the thermal lightcurve exercise 
as we did for the PUI data. The predicted 70.0 fim lightcurve am- 
plitude is about 1 mJy peak-to-peak, while the flux error of the 
best 70.0 yum data point is already 1 .5 mJy. 

4.5. Size 

The unknowns in surface roughness influenced the optimum 
thermal inertia slightly (see above), but the size and albedo val- 
ues did not change significantly. We repeated the exercise shown 
in Fig.[3]with a the low surface roughness (rm.s. -slope p = 0.2), 
which produced the best match to the Spitzer-IRS data using a 
spherical and the ellipsoidal shape model specified above. The 

values are now a few percent lower due to the better fit to the 
shortest wavelength data and, in the case of the ellipsoidal shape 
model, also a better match to the PUI data. The ;t'^-minimum 
is pointing towards a thermal inertia of r=650 Jm"^s"°^K"\ re- 
sulting in an equivalent diameter of an equal volume sphere of 
D,// = 494+j|m when we consider the full 3-cr confidence level 
for the thermal inertia. The VISIR data are better matched by a 
relatively high roughness. The derived corresponding values are 
then D,ff = 508!^; m (atp = 1.0, Fop,,,™,,, = 850 Jm-^s-O-SK-'). 
The equivalent size of 1999 RQ36 (of an equal volume sphere) 
is therefore in the range between 480 to 511m. The true size 
error is probably a bit larger because the equator-on observing 
geometry does not allow us to "see" the pole regions in the mid- 
to-far-IR thermal emission. 

The derived diameter is significantly smaller than the 610 m 
derived from Spitzer measurements alone (Emery et al. 2010). 
Such a large size is clearly not compatible with the full set 
of thermal measurements. Emery (private communication) con- 
firms that the published size was produced by a simple thermal 
model (NEATM, Harris I1998I I. which is known to overestimate 
the diameters of asteroids by 15-20% especially when the ther- 
mal dataset is limited to only large phase angles (Harris 1213061) . 

It is interesting to note that Hudson et al. (I2000I I used the 
first Goldstone radar images -taken in 1999- to derive a diame- 
ter estimate of ~500m, very close to our value. Their prelimi- 
nary results indicate a remarkably featureless surface (down to 
the resolution limit of 19 m) and an almost spherical object, but 
they derived a wrong spin period of ~2 hr. Later on, Nolan et al. 
dSOOTl l used multiple radar data from 1999 and 2005, some of 
them taken with unusually good aspect coverage, and conclude 
from the combined dataset that 1999 RQ36 must be an irregular 
spheroid with about 580 m in diameter, more than 15% larger 
than our value. Only the recent work by Nolan et al. (I2012I I has 
brought the size back to the 500 m value. 

4.6. Geometric albedo 

The derived albedo is in the mid-range of values suggested for 
B-type asteroids (Tholen & Barucci 1989). But, as noted be- 
fore, the albedo value depends directly on the H-magnitude. 
Our solution of 0.045;^q q[2 already includes the error of +0.3 
in H-magnitude. If it turns out that the opposition effect of 
1999 RQ36 is steeper or shallower than assumed for calculating 
of the H-magnitude, then the derived geometric albedo would 
also change. Hergenrother et al. (120121 ) provide, in addition to 
the above H-magnitude, Hv and G by strictly following the lAU 



H-G photometric system for airless bodies (Bowell et al. 119891) 
and find Hi/=19.90+0.10mag and G=-0.14+0.02. The corre- 
sponding radiometric size would not change, but the geometric 
albedo could then be as high as pv=0.085. The low albedos for 
near-Earth objects are believed to be connected to a primitive 
volatile-rich surface composition (Fernandez et al. 120051) . But a 
dedicated orbital evolution and thermal history study by Delbo 
& Michel (120111) has shown that 1999 RQ36 has approached 
the Sun in the past to regions of q < 0.8 (corresponding to 
about 400 K at the sub-solar point) with an 80% probability and 
with a lower probability even to q < 0.6 (480 K at the sub-solar 
point). Such high temperatures alter the pristine properties of the 
surface material, because some of the expected primitive com- 
pounds already break up at moderate temperatures above 300 K, 
and the surface volatiles might be completely lost (e.g., Marchi 
et al. 120091 ). But the temperature maxima at 3-5 cm below the 
surface are already lower by 100 K, and organics have very likely 
been protected from thermal break-up (Delbo & Michel [201 11) . 

4.7. TPM concept 

The calculations for the TPM beaming function for a rough sur- 
face are, strictly speaking, only applicable to the sunlit part of 
the asteroid (Lagerros 1998). A change in roughness changes 
the radiation field. A rougher surface has the tendency to emit 
or "beam" more in the solar direction at the expense of emission 
into larger phase angles. And this calculation is only done for 
the sunlit part of the surface. But even at large phase angles, the 
observed mid- and far-IR thermal emission is still dominated by 
the emission from the sunlit part of the object (^T"*), the emis- 
sion from the dark and cold part beyond the terminator is negli- 
gible with or without roughness. Only for sub-millimetre/milli- 
metre/centimetre wavelengths observations at very large phase 
angles above 90° are the night side contributions not negligible 
anymore, and these assumptions might turn out to be problem- 
atic and might result in wrong predictions. Muller et al. (20051) 
used observations up to 110° phase angle to derive highly accu- 
rate size and albedo values for 25143 Itokawa. The correspond- 
ing roughness and thermal inertia settings (similar to the values 
found here) explain the observations taken at 27° and 54° with 
a similar accuracy to the observations taken at 108-110°, giv- 
ing us confidence that our TPM setup is appropriate also for the 
1999 RQ36 observations. 

5. Conclusions 

The Herschel DDT measurements of 1999 RQ36, combined 
with ESO- VISIR DDT measurements and data from the Spitzer 
Space, have allowed us to conclude that this asteroid has a diam- 
eter of an equal volume sphere in the range 480 to 5 1 1 m, consid- 
erably smaller than previously estimated radar and radiometric 
values (Emery et al. 120101 Nolan et al. 2007), but in good agree- 
ment with the recent mean diameter given by Nolan et al. ( 120121) 
of 493 +20 m. The thermal lightcurve variations are best-fit by 
an ellipsoidal shape model with a best fit a/b-ratio of of 1.04, 
corresponding to an ellipsoidal body with 509 x 489 x 489 m (+ 
about 10 m in a- and b- and an unknown error in c-dimension). 
1999 RQ36 has a geometric albedo pv - 0.045^^^-^^\l, in the nom- 
inal range for primitive volatile-rich NEAs of this type, and has 
a thermal inertia r=650 + 300 Jm"^s"°-^K"', similar to Itokawa, 
a rubble pile asteroid. Our results further serve to confirm that 
the asteroid has a retrograde rotation with a spin vector between 
-70 to -90°. The inconsistency between the flux value derived 
from VISIR observations and the Spitzer-IRS dataset, taken at 
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very different observing geometries, suggests there may be large 
variations in roughness on the surface. If the motion of surface 
material or regolith turnover processes are responsible for differ- 
ent terrains, then one can expect to find regions on the surface 
with very fresh material, at least at subsurface layers at several 
centimetres depth. The results of this work on the hypothesis 
that two different types of terrains exist with different levels of 
roughness and the possibility of having part of the surface with 
fine dust and pebbles is very important for the OSIRIS-REx mis- 
sion sample mechanism strategy. A global mapping during the 
presampling phase will then allow the OSIRIS-REx team to se- 
lect the most optimal area for collecting pristine and organic -rich 
material. 

With this project we explore possibilities for thermophysi- 
cal model techniques for an individual target that will soon be 
seen and characterised in detail. It will allow us to verify our 
findings on physical and thermal properties, as well as the pro- 
posed different surface terrains that are needed to explain the 
derived roughness and thermal inertia. Radiometric techniques 
are very powerful for deriving highly reliable size and albedo 
information, but signatures of surface texture, shape, spin-state, 
and thermal properties are also included in the object's thermal 
emission. In fact, because thermal data is (or will soon be) avail- 
able for many thousands of asteroids (IRAS, MSX, ISO, Akari, 
Spitzer, WISE, ground-based mid-IR/sub mm/mm programmes, 
etc.), the results from 1999 RQ36 can be considered to be the key 
to allowing us to transfer our model techniques to many other 
targets that will not be visited by spacecraft in the near future, 
but for which similar important questions can now be answered. 

Acknowledgements. HSpot & HIPE are joint developments by the Herschel 
Science Ground Segment Consortium, consisting of ESA.the NASA Herschel 
Science Center, and the HIFI, PACS, and SPIRE consortia. We gratefully ac- 
knowledge the awarding of director discretionary time from the Herschel Space 
Observatory and ESO to support the work published here. We would also like to 
acknowledge the generous support by Mario van den Ancker (ESO) for instal- 
lating and running of the VISIR pipeline and the support by Eric Pantin (CEA) 
and Ralf Siebenmorgen (ESO) for VISIR calibration and eiTor calculation is- 
sues. Support for the Spitzer data reduction was provided by Sonia Fornasier, 
Josh Emery, and Jeffrey Van Cleve. It is also our pleasure to thank the Herschel 
Science Centre Mission Planning Team for their efforts in the optimal scheduling 
of these observations. Cs. K. and A. P. are supported by the ESA Grant PECS- 
98073 and by the Bolyai Research Fellowship of the Hungarian Acedemy of 
Sciences. We would also like to thank the anonymous referee for a very thor- 
ough and very constructive review of the manuscript. 



Krugly, Y. N.. Belskaya. I. N., Shevchenko, V. G. et al. 2002. Icarus, 158, 293- 
304 

Lagage, R O., Pel, J. W., Authier, M. et al. 2004, The Messenger 1 17, 12-16 
Lagerros, J.S.V, 1996, A&A, 310, I0II-I020 
Lagen-os, J.S.V, 1997, A&A, 325, 1226-1236 
Lagen-os, J.S.V, 1998, A&A, 332, 1123-1132 

Lauretta, D. S., Drake, M. J., Binzel, R. P, et al. 2010, M&PSA, 73, 5153 
Lauretta, D. S. & The OSIRIS-REx Team 2012, LPS XXXXIII, Abstract #2491 
Marchi, S., Delbo, M., Morbidelli, A. et al. 2009, MNRAS 400, 147-153 
Delbo, M. & Michel, P 201 1, ApJL, 728, 42 

Milani, A., Chesley, S. R., Sansaturio, M. E. et al. 2009, Icarus 203, 460 

Mommert, M., Harris, A. W., Kiss, C. et al. 2012, A&A, accepted 

Muller, T. G. & Lagerros, J. S. V. 1998, A&A, 338, 340-352 

Muller, T. G., Lagerros, J. S. V, Burgdorf, M. et al. 1999, in The Universe as 

Seen by ISO. Eds. R Cox & M. F Kessler, ESA-SP 427, 141-144 
Muller, T G. & Lagerros, J. S. V. 2002, A&A, 381, 324-339 
Muller, T. G. 2002, M&PS, 37, 1919 

Muller, T. G., Sekiguchi, T, Kaasalainen, M. et al. 2005, A&A, 443, 347-355 
Muller, T. G., Durech, J., Hasegawa, S. et al., 201 1, A&A, 525, A145 
Nolan, M. C, Magri, C, Ostro, S. J. et al. 2007, AAS DPS # 39, #13.06, BAAS 
39, 433 

Nolan, M. C, Magri, C, Benner, L. A. M. 
et al. 2012, ACM 2012 abstract # 6 345, 
http: //www. Ipi .usra.edu/meetings/acra2812/pd£/6345 .pdf 

O'Rourke, L.. Muller, T. G., Valtchanov, I. et al. 2012, P&SS, accepted 

Ostro, S. J., Margot, J.-L., Benner, L. A. M. et al. 2006, Science, 314, 1276-1280 

Ott, S. 2009, Herschel Interactive Processing Environment, HIPE, Proceedings 

of the Astronomical Data Analysis Software and Systems XIX Conference, 

434. 139-142 

PACS photometerPoint Source Flux Cahbration, Muller et al. 2011, PICC-ME- 
TN-037, Version 1.0, retrieved 29/Feb/2012; technical report available via 



URL http://herschel.esac.esa.inl/twiki/bin/view/Public/PacsCalibrationWeb/pacs_bolc 
PACS Photometer Passbands and Colour Conection Factors for Various 
Source SEDs, Muller, T, Okumura, K., Klaas, U. 2011, PICC-ME- 
TN-038, Version 1.0, retrieved 29/Feb/2012; technical report available via 



URL http://herschel.esac. esa.int/twiki/pub/Public/PacsCalibrationWeb/cc_report_vl.p 
Pilbratt, G. L., Riedinger, J. R., Passvogel, T. et al. 2010, A&A, 518, LI 
Poglitsch, A., Waelkens, C, Geis, N. et al. 2010, A&A, 518, L2 
Santos-Sanz, P., Lellouch, E., Fornasier, S. et al. 2012, A&A, accepted 
Schiitz, O. & Sterzik, M. 2005, in High Resolution Infrared Spectroscopy in 

Astronomy, Proceedings of an ESO Workshop held at Garching, Germany, 

18-21 November 2003. Kaufl, Siebenmorgen, Moorwood (Eds.), 104 
Tholen, D. J. & Baracci, A. M. 1989, Asteroids II; Proceedings of the 

Conference, Tucson, AZ, Mar. 8-11, 1988 (A90-27001 10-91), University 

of Arizona Press, 1989, p. 298-315 
Vilenius, E., Kiss, C, Mommert, M. et al. 2012, A&A, submitted 
Vokrouhlicky, D., Nesvorny, D., Bottke, W. F 2003, Nature, 425, 147 
Walsh, K. J., Richardson, D. C. & Michel, P 2008, Nature 454, 188-191 
Werner, M. W., Roelhg, T. L., Low, F J. et al. 2004, ApJS, 154, 1 



References 

Binzel, R.P, Morbidelli, A., Merouane, S. et al. 2010, Nattjre, 463, 331 
Bowell, E. 1989. Asteroids II; Proceedings of the Conference, Tucson, AZ, Mar. 

8-II, 1988 (A90-27001 10-91), University of Arizona Press, 1989, p. 524- 

556 

Campins, H., Morbidelh, A., Tsiganis, K. et al. 2010, ApJ, 721, 53 
Cellino, A., Zappala, V., Doressoundiram, A. et al. 2001, Icarus, 152, 225-237 
Clark, B. E., Binzel, R. R, Howell, E., et al. 201 1, Icarus, 216, 462-475 
Delbo, M., Dell'Oro, A., Harris, A. W. et al. 2007, Icarus, 190, 236 
Dombard, A. J., Bamouin, O. S., Prockter, L. M. & Thomas, P C. 2010, Icarus, 
210, 713 

Emery, J. P, Fernandez, Y. R., Kelley, M. S., et al. 2010, 41st LPSC, LPI contri- 
bution No. 1533, 2282 
Fernandez, Y. R., Jewitt, D. C, Sheppard, S. S. 2005, AJ, 130, 308 
HaiTis, A. W. 1998, Icarus 131, 291 

Harris, A. W. & Lagerros. J. S. V. 2002, Asteroids III, W. F Bottke Jr., A. Cellino, 
P. Paolicchi, and R. P. Binzel (eds). University of Arizona Press, Tucson, 
p.205-218 

Harris, A. W. 2006, in Asteroids, Comets, Meteors, eds. D. Lazzaro, S. Ferraz- 

Mello & J. A. Fernandez, Proc. lAU Symp., 229, 2005 
Hergem-other, C. W., Nolan, M., d'Aubigny, C. et al. 2012, LPS XXXXIII, 

Abstract #2219 

Hudson, R. S., Ostro, S. J., Benner, L. A. M. et al. 2000, AAS DPS # 32, #07.10, 
BAAS 32, 1001 



9 



